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I.  GRITS  RAL 


An  objective  analysis  using  normal  mode  functions  evaluated  at 
irregularly  spaced  locations  was  designed  and  developed  to  provide  initial 
fields  for  the  AFGL  global  spectral  model.  The  software  for  full  forecast 
and  analysis  cycles  was  implemented  on  the  Air  Force  Weapons  Laboratory 
(AFWL)  Cray-1  system.  The  objective  of  this  user  manual  is  to  provide  an 
overview  of  the  flow,  input  requirements,  and  output  capabilities  of  the 
objective  analysis  scheme.  Technical  details  are  discussed  in  Halberstam 
et  al.,  1984. (1> 


1.  Halberstam,  M. ,  C.  Johnson,  D.  C.  Morquist,  S.-L.  Tung,  1984:  Two 
Methods  of  Global  Data  Assimilation.  AFGL-TR-84-0260,  Contract  F19628-82- 
C-0023,  Systems  and  Applied  Sciences  Corporation. 
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SYSTEM  FLOWCHART 


System  procedures  for  each  analysis  cycle  without  iteration: 


System  procedures  with  three  iterations  for  each  analysis  cycle: 


III.  DOCUMENTATION  OF  SEQUENTIAL  MAIN  PROGRAMS 


All  the  programs  developed  for  the  analysis  cycle  are  documented  in 
this  section.  This  documentation  is  written  in  the  following  form  for  each 
program: 

PROGRAM:  (Program  name) 

ABSTRACT :  (Abstract  of  the  program) 

MAIM  SUBROUTINE:  (Name  of  main  subroutine  or  calling  subroutine 

of  the  program) 

SUBROUTINES  CALLED:  (Name  of  subroutine  used  in  the  program) 

INPUT  FILES:  (Structure  and  variables  for  input  files  in  this 
program) 

OUTPUT  FILES:  (Structure  and  variables  for  output  files  in  the 
program) 

DESCRIPTION  OF  MAIN  SUBROUTINE:  (Detail  of  the  coding) 


PROGRAM; 


C2POSTPTZ 

C11POSTPTZ 


ABSTRACT:  These  programs  are  set  up  to  process  the 

spectral  data  from  the  forecast  model  to 
form  a  first  guess  grid  field  for  input  to 
program  C8GRIDF  to  reconstruct  a  new 
analyzed  field.  The  temperature  field  is 
converted  to  height  field  by  using  Sela's 
method. 


MAIM  SUBROUT IHE :  FIELDC(C2POSTPTZ)- -input  formatted  spectral 

data  (coded  form) 

FI ELDB (Cl 1P0STPTZ ) -  input  unformatted 

spectral  data 
(binary  form) 

SUBROUTIMES  CALLED:  LEGSUM 

FFT1 

PMNS 

BSCST 

SPTOGP 

UMVM 

SIMMC 


OUTPUT  FILES: 


1.  FORTRAN  UNIT  2  -  First  guess  grid  field 

Record  #  Description 

1  to  (NLAT*NHEM)  WTR(MP)  =  Surface  pressure 

(NLAT*NHBM+1)  to  (KP*3+1*NLAT*NHEM)  WTR(MP)  =  Written  in  sets  of  3  (U, 
NLAT=Total  number  of  latitudes  V,  and  heights)  for  each  sigma 
NHEM=Total  number  of  hemispheres  layer 

KP=Total  number  of  sigma  layers  MP  =  Total  number  of  longitudes 

DESCRIPTION  OF  MAIN  SUBROUTINES: 

Designators  Text 

CA  Call  BSCST  to  set  up  constants,  resolutions 

and  sigma  structure. 

Call  CFFTI  to  set  up  Fourier  transformation, 
the  dimension  of  WSAVE  is  MP*4+15. 

CB  Read  in  the  12  hours  forecast  spectral  data. 

CC  Calculate  geopotential  by  using  Sela's 

hydrostatic  equations. 

CD  Start  Loop  (DO  6000)  over  number  of 

latitudes  to  transform  dependent  variables 
from  spectral  (spherical  harmonic)  space  to 
grid  space. 

Call  SIGNMC  to  reset  the  data  from  AFGL 
sigma  structure  (K=l  for  top),  to  NMC  sigma 
structure  (K=l  for  bottom). 

CE  Reset  the  data  structure  for  output  file 

(tape  2) . 


PROGRAM: 


C4RFMPRO 


ABSTRACT:  This  program  is  designed  to  reform  the  residuals  file  and 

observations  file  from  the  output  of  program  C3ASAPRES. 

Both  temperature  residuals  and  temperature  observation 
profiles  are  converted  to  height  residuals  and  height 
profiles  by  using  the  method  of  Sela  (1982). (2) 

Zero  residuals  are  filled  in  Where  data  are  missing.  The 
sounding  is  rejected  if  the  data  from  all  layers  are 
missing.  The  update  surface  pressures  are  computed  by 
using  a  quadratic  relationship.^^  Four  corners  method 
is  used  to  insert  the  first  guess  for  data-void  area.^1* 

MAIM  SUBROUTINE:  RFMZERO 

SUBROUTINES  CALLED:  DTTODZ 

GETZERO 

TTOZ 

GBTPS 

SETSIG 


INPUT  FILES: 

1.  FORTRAN  UNIT  1  -  Buddy-checked  residuals  and  observations  data  at 

observation  sites  (FGGE  II)  which  output  from 
C3ASAPRES  program 

Records  »  Description 

1  NOBS  *  Number  of  soundings 


Even  of  (2  to  NOBS  +  1)  RLAT.RLON, JID, ZSTAR, (UBD( J) , J=1 ,97) 

RLAT  =  Latitude  of  the  observation  site 
RLON  =  Longitude  of  the  observation  site 
JID  =  Soundings  ID 
ZSTAR  =  Surface  height 
0BD(1)  =  Surface  pressure 
0BD(2)  =  Height  residual  at  level  2 
OBD ( 3 )  =  Quality  flag  for  height 
0BD(4)  =  U  wind  residual  at  layer  1 
OBD ( 5 )  =  V  wind  residual  at  layer  1 
0BD(6 )  =  Quality  flag  for  wind 
0BD(7)  =  Moisture  residual  at  layer  1 
0BD(8)  =  Quality  flag  for  moisture 
0BD(9)  =  Temperature  residual  at  layer  1 


2.  Sela,  J.,  1982:  The  MMC  Spectral  Model.  NOAA  Technical  Report  NWS-30, 
National  Meteorological  Center,  Washington,  D.C. 


Records  # 


Odd  of  (2  to  NOBS+1) 


2.  FORTRAN  UNIT  10  - 
equations 

Record  # 

1 

2 

OUTPUT  FILES: 

1.  FORTRAN  UNIT  5  - 

Record  # 

(K-1)*KP+1 

K  =  Kth  sigma  layer 
KP  =  Total  number  of 
sigma  layers 

(K-l)*KP+2 

(K-l)*KP+3 

(K-l)*KP+4 


Description 


LL=8*(K-1)+1 

OBDCLL+1)  =  Height  residual  at  level  K+l 
0BD(LL+2)  =  Quality  flag  for  height 
0BD(LL+3)  =  U  residual  at  layer  K 
0BD(LL+4)  =  V  residual  at  layer  K 
0BD(LL+5)  =  Quality  flag  for  wind 
0BD(LL+6 )  =  Moisture  residual  at  layer  K 
0BD(LL+7)  =  Quality  flag  for  moisture 
0BD(LL+8)  =  Temperature  residual  at  layer  K 


(TT( J) , J=1 ,KP) , Xl( J) , J=1 ,KP) , 
<X2(J),J=1,KP),<X3(J),J=1,KP), 

(TFG(J) , J=1,KP) 

TT  =  Temperature  observation  profile 

XI  =  U  wind  observation  profile 

X2  =  V  wind  observation  profile 

X3  =  Specific  humidity  observations  profile 

TFG  =  First  guess  temperature  profile 


matrix  and  C  matrix  for  Sela's  hydrostatic 


Description 


<AH<  I ,  J) ,  I»1 ,  KP) ,  J=1 ,  JCP> 

AM  i  A  matrix 
(CM(I,J),I«1,KP),J=1,KP) 

CM  =  C  matrix,  (the  inverse  of  AM)*  B  matrix 


Reorganized  residual  data  of  U,V  wind  and 
composite  variable  P  with  equal  number  of  data 
points  at  each  layer 

Description 

N  =  Number  of  soundings  for  U  residuals 


(OBSLAT(I) , I=1,N) 

Latitudes  for  U  residuals 

(OBSLONC I ) , 1=1 ,N) 

Longitude  for  U  residuals 

(X(I),I=1,N) 

U  wind  residuals 
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(K-l)*KP+5 


H  =  Humber  of  sounding  for  V  residuals 


(K-l)*KP+6 

(K-1)*KP+? 

(K-l)*KP+8 

(K-l)*KP+9 

(K-1)*KP+10 

(K-1)*KP+11 

(K-1)*KP+12 

FORTRAN  UNIT 

Record  # 

(K-1)*KP+1 

(K-l)*KP+2 

(K-l)*KP+3 

CK-1)*KP+A 

(K-l)*KP+5 

(K-l)*KP+6 

(K-l)*KP+7 

(K-l)*KP+8 


(OBSLAT(l) , 1=1 ,N) 

Latitude  for  U  residuals 

(OBSLON(I) ,1=1 ,M) 

Longitude  for  V  residuals 

(X(I)  ,1=1  *N) 

V  wind  residuals 

H  =  Humber  of  soundings  for  P 
(composite  variable  residuals) 

(OBSLAT(I) , 1=1 ,N) 

Latitude  for  P  residuals 

(OBSLOH(I) , 1=1 ,N) 

Longitude  for  P  residuals 

(X(I) ,1=1, H) 

P  residuals 


15  -  Reorganized  observation  data  of  U,V  wind  and 

temperature  with  equal  number  of  data  points 

at  each  layer 

Description 

N  =  Number  of  soundings  for  U 

(OBSLAT(I) ,1=1, N) 

Latitudes  for  U  observations 

(OBSLON(I) ,1=1 |H) 

Longitudes  for  U  observations 

(X(I),I=1,N) 

U  wind  observations 

H  =  Number  of  soundings  for  V  observations 

(0BSLAT(2) , 1=1 ,N) 

Latitudes  for  U  observations 

(OBSLON(I) ,1=1 ,N) 

Longitudes  for  V  observations 

(X(I) ,1=1 ,N) 

V  wind  observations 
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(k-1)*kp+9 


II  =  Number  of  soundings  for  temperatures 


(K-1)*KP+10 

(K-1)*KP+11 

(K-1)*KP+12 


(OBSLAT(I) ,1=1  ,N) 

Latitudes  for  temperatures 

(OBSLON(I) , 1=1 ,N) 

Longitude  for  temperatures 

(X(I) ,1=1 ,N) 


Temperature  observations 
DESCRIPTION  OF  MAIN  SUBROUTINE: 


Designators  Text 

CA  Read  in  total  number  of  soundings. 

CB  Set  up  sigma  structure. 

(K=l  for  bottom) 

CC  Read  in  matrix  A  and  matrix  C. 

CD  Input  residuals,  surface  height,  observations  and 

first  guess  temperature  on  each  sigma  layer. 

CE  Check  the  missing  data.  If  the  data  of  all  the 

layers  are  missing,  the  sounding  will  be  dropped. 

If  the  data  are  found  on  at  least  one  layer,  the 
sounding  will  be  reformatted  and  stored.  The 
missing  residuals  will  be  replaced  by  zeros  and 
the  missing  observation  temperatures  will  be 
replaced  by  the  first  guess  temperatures . 

CF  Sela's  method  is  used  in  both  subroutines  DTTODZ 

and  TTOZ  to  convert  temperature  residuals  te 
height  residuals  and  temperature  profiles  to 
height  profiles.  Subroutine  GSTPS  is  called  to 
compute  the  surface  pressure  8t  observation  sites 

CG  The  residuals  of  composite  variable  are  confuted 

from  the  update  surface  pressure  residuals  and 
height  residuals. 

AP  =  AZ*9 . 8+R*To*AQ 
R  =  287.05 

To  =  300°K  basic  state  temperature. 

CH  The  useful  data  are  stored  in  the  scratch  files 

tape  2  (residuals)  and  tape  12  (observations). 
Subroutine  GETZERO  is  called  to  apply  four 
corners  method  to  fill  in  first  guess  data  for 
data-void  area. 


PROGRAM: 


C5VERMODES 


ABSTRACT:  This  program  is  designed  to  read  data  for 

the  sigma  layers,  then  project  the  data  on  desired 
vertical  modes. 


MAIM  SUBROUTINE:  VERHOD 


SUBROUTINES  CALLED: 

VERT 

INPUT  PILES: 

1.  FORTRAN  UNIT  1  - 

U,  V,  P  residuals  at  sigma  layers 

(See  output  file  of  program  C4RFMPRO) 

2.  FORTRAN  UNIT  5  - 

Vertical  modes 

Record  # 

Description 

1 

((ErGG(I,J),I=l,KLBV),J=l, MODES), 

( (EIGT(I , J) ,1=1 ,KLEV) , J=1 .MODES) , 

(GH(I) ,1=1 .MODES) , (DOTPRO(I) ,1=1 .MODES) , 

To(I) ,I=1,KLEV) 

KLEV  =  Number  of  sigma  layers 

MODES  =  Number  of  vertical  modes 

BIGG  =  Eigenvectors  of  matrix  G*3* 

BIGT  =  Eigenvectors  of  matrix  GT 

GH  =  Negative  of  the  real  eigenvalues  of 
matrix  G 

DOTPRO  =  Reciprocal  of  the  inner  product 
of  eigenvectors  of  G  and  GT 

To  =  A  vector  of  basic  state  temperature 

OUTPUT  FILES: 

1.  FORTRAN  UNIT  10 

-  Data  projections  on  vertical  modes 

Record  # 

Description 

<IM0D-1)*4+1 

KOBS , IMODES , I MOD 

KOBS  =  Number  of  soundings 

IMODES  =  Number  of  vertical  modes  applied 

IMOD  =  Mode  index  (1  to  IMODES) 

(IMOD-l)*4+2 

(RANLAT( J) , J=1 ,K0BS) 

Latitudes  for  soundings 

3.  Ballish,  A.  B. , 
NMC  Spectral  Model. 

1980:  Initialization.  Theory  and  Application  to  the 
Ph.  D.  Thesis,  University  of  Maryland,  151  pp. 
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Record  » 

(  MOD- 1 )  *4+3 

(  MOD- 1)  *4+4 


Description 

(RAHLON(J) , J=1 , KOBS) 

Longitude  for  soundings 

(OBS(J,IMOD) , J=KOBS) 

Data  projections  of  mode  MOD 


DESCRIPTION  OF  MAIM  SUBROUTINE: 

Designators  TEXT 

CA  IITPUT  vertical  modes  from  Tape  5. 

INPUT  residuals  from  Tape  1. 

CB  Call  subroutine  VERT  to  do  vertical  projection. 

The  projection  is  located  at  observations  site 
for  each  MOD,  up  to  MODES. 

CC  Write  the  projections  to  Tape  10  for  output. 

Print  out  some  values  for  checking. 


PROGRAM:  C6DIXBXPNO 

ABSTRACT:  This  program  is  designed  to  perform  the 

objective  analysis  on  irregularly  spaced 
observation  points  by  fitting  normal  mode 
functions  to  the  residuals  (observations  minus 
first  guess  values)  with  "finality" 
procedure. (D  The  residuals  are  first 
projected  onto  vertical  modes  and  a  set  of 
horizontal  normal  functions  for  each  vertical 
mode  is  evaluated  at  the  same  geographic 
locations  as  the  observations.  A  set  of 
coefficients  is  determined  by  fitting  the 
vectors. 

MAIM  SUBROUTINE:  PCOMP 

SUBROUTIMBS  CALLED:  GETTH 

GETOBS 

PUNCH 

IHPUT  FILES: 

1.  FORTRAN  UNIT  1  -  DATA  projections  on  vertical  modes 

(See  output  files  of  program  C5VERMODES) 

2.  FORTRAN  UNIT  2  -  Normal  mode  functions  for  every  one  degree  latitude 

interval 


Description 

GH  s  Geopotential  of  vertical  mode  IMOD 


IS , IWG , LR , LALL , IAS 


IS  a  Zonal  wave  number 
IWG  s  Gravity  wave  number 
LR  3  Rossby  wave  number 
LALL  =  Accumulated  wave  number 
IAS  3  index  for  symmetric  (=1)  or  antisymmetric 
(=2) 


(FH(K, 1 ) ,K=1 ,NLAT) 

Normal  mode  functions  for  variable  U  of  zonal 
wavenumber  IS  and  frequency  index  LALL  at  every 
1  degree  latitude 

(FH(Kt2) ,K=1 ,NLAT) 

Normal  mode  functions  for  variable  V  of  zonal 
wavenumber  IS  and  frequency  index  LALL  at  every 
1  degree  latitude 
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Record 


10  to  IM0DES*NS*NLRM+1 
NS  =  Total  number  of 
zonal  waves 


Description 

(PH(K,3),K*1,HLAT) 

Normal  mode  functions  for  variable  P  of  zonal 
wavenumber  IS  and  frequency  index  LALL  at  every 
1  degree  latitude 

IS,IWG,LR,LALL,IAS 

(FH(K,1) , K=1,NLAT) 

(FH(K,2) , K=1,NLAT) 

(FH(K,3),K=1,NLAT) 

Pattern  of  record  #2  to  record  #5  is  repeated 
for  all  the  zonal  wavenumbers  and  all  the 
frequency  indices  of  each  vertical  mode  IMOD. 
IMODES=8  is  used  for  case  of  12  sigma  layers 
model. 


3.  FORTRAN  UNIT  3  -  Total  number  of  gravity  waves  and  Roseby  waves 

to  be  used  of  each  zonal  wave  at  each  vertical  mode 


Record 


1  to  IM0DES*NS 


Description 

NLR  =  Number  of  frequency  indices  for  each 
zonal  wavenumber  IS  and  each  vertical  mode  IMOD 
NL  =  Maximum  of  NLR 


OUTPUT  FILES: 

1.  FORTRAN  UNIT  20  -  Resulting  coefficients  from  analysis 


Record 


1  to  T MODES 


Description 


A(NL,NS,2)  3  Coefficients  for  each  vertical  mode 


DESCRIPTION  OF  MAIN  SUBROUTINE: 


Designators 


Input  resolution  and  constants. 

In  loop  (DO  1000),  vertical  modes  and  normal 
mode  function  are  read  in  from  Tape  2. 

For  each  vertical  mode,  subroutine  GETOBS  is 
called  to  get  data  projections  of  this  vertical 
mode  and  subroutine  GETTH  is  called  to  evaluate 
the  normal  mode  functions  on  observation  sites. 

"Finality"  procedure  is  performed  in  inner  loop 
(DO  5)  on  all  the  observation  points  for  each 
zonal  wavenumber  (DO  8)  and  each  frequency  index 
(DO  7). 
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CD  The  coefficient  array  AdtL.MS^)  is  written  to 

Tape  20. 

CB  The  differences  between  the  input  data 

projections  and  the  computed  values  from 
coefficients  at  the  observation  sites  are 
printed.  The  RMS  errors  are  computed  for  each 
variable  (U,  V,  P)  at  each  vertical  mode. 


PROGRAM: 


C7DIXARL 


ABSTRACT:  This  program  is  designed  to  use  the  coefficients 

output  from  C6DIXEXPM0  to  determine  the  values 
of  variables  at  regular  grid  points (2.5*  X 
2.5*  global  grid). 

MAIM  SUBROUTINE:  DIXAHL 

SUBROUTIME  CALLED:  Hone 

IMPUT  FILES: 

1.  FORTRAM  UNIT  3  -  Total  number  of  gravity  waves  and  Rossby  waves  to  be 

used  for  each  zonal  wave  at  each  vertical  mode 
(See  FORTRAN  UNIT  3  input  file  for  program 
C6DIXEXPN0 . ) 

2.  FORTRAN  UNIT  4  -  Normal  modes  functions  for  every  2.5  degree  latitude 

interval 

(The  file  structure  is  the  same  as  the  structure  of 
input  file  FORTRAN  UNIT  2  for  program  C6DIXEXPMO.) 

3.  FORTRAN  UNIT  5  -  Vertical  modes 

(See  input  file  FORTRAN  UNIT  5  for  program 
C5VERM0DES .) 

4.  FORTRAN  UNIT  20  -  Analyzed  coefficients  from  the  output  of  program 

C6DIXBXPM0 

(See  FORTRAN  UNIT  20  of  program  C6DIXEXPNO.) 


OUTPUT  FILES: 

1.  FORTRAN  UNIT  30  -  Analyzed  residuals  for  each  variable  (U,  V,  P)  at 

regular  global  grids  (2.5*  X  2.5*)  on  sigma  layers 


Record 


Description 


1  to  KP*MFUM*LHAF  ( (X(I , J, ILAT , IFUN) , 1=1 ,HL0N) , J=1 , 2) 

This  record  structure  is  repeated  for  LHAF  latitudes 
and  NFUN  variables  and  KP  sigma  layers 
NLON  =  Number  of  longitude  points 

(NL0N=144  for  every  2.5  degree) 

J  =  Hemisphere  indicator 

(J=l  for  northern  hemisphere, 

J=2  for  southern  hemisphere) 


Record  #  Description 

ILAT  =  Latitude  index 

(ILAT  =  1  to  LHAF) 

I FUN  =  Variable  index 

(IFUMsl  for  U  wind  residuals, 

IFUN»2  for  V  wind  residuals, 

IFUN=3  for  surface  pressure  residuals  ) 

LHAF  =  Number  of  latitudes  for  one  hemisphere 
(LHAF- 3 7  for  every  2.5  degree  latitude) 

NFUN  =  Number  of  variables  for  processing  (NFUN=3), 
i.e.,  U,  V,  P  (composite  variable) 
or  U,  V,  Q  (surface  pressure) 

DESCRIPTION  OF  MAIN  SUBROUTINE: 

Designators  Text 

CA  Input  for  resolutions  and  constants. 

CB  Set  up  sigma  structure  and  basic  state  temperature. 

CC  Initialize  the  variable  array  and  input  the  vertical 

modes  from  Tape  5. 

CD  Compute  the  scalar  product  of  the  vector  function  of 

the  sigma  coordinates  with  the  vertical  mode 
eigenvector. 

CE  Get  coefficients  SGAM(NL,NS,2)  from  Tape  20  and  get 

normal  mode  functions  from  Tape  4. 

Residuals  for  variables  (U,  V,  P)  on  a  2.5*  X  2.5° 
global  grid  for  each  vertical  mode  IHOD  are 
determined  (GRID)  and  stored  in  scratch  file  Tape  1. 

I MOD  ranges  from  1  to  INODES  (IM0DBS=8  for  KP=12 
sigma  layers). 

CF  At  each  grid,  values  on  the  vertical  modes  are 

converted  to  the  values  on  the  sigma  layers  and 
composite  variables  are  converted  to  heights. (3) 
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PROGRAM: 


C8GRIDP 


ABSTRACT:  This  program  is  set  up  to  reconstruct  the  grid  field 

from  analyzed  residuals  field  (2.5°  X  2.5°).  The 
updated  surface  pressure  field  is  computed  from 
updated  height  field  by  using  a  quadratic 
relationship.  The  height  field  is  converted  to 
temperature  field  by  using  the  Flattery 
method. ID  All  the  variables  are  interpolated  to 
updated  sigma  layers. 


MAIN  SUBROUTINE: 

GRIDF 

SUBROUTINES  CALLED: 

PSZTOT 

FMTTRA 

SETSIG 

TOSIG 

FLATZT 

IMINV 

INPUT  FILES: 

1.  FORTRAN  UNIT  1  - 

First  guess  field  (2.5°  X  2.5°) 

Record  # 

Description 

1  to  74 

Q(LON)  =  Surface  pressure 

LON  a  Ntunber  of  points  for  each  latitude  circle 
(L0N*144  for  every  2.5*  longitude) 

ODD  record  #  is  for  northern  hemisphere  and  even 
record  #  is  for  southern  hemisphere,  running  from 
equator  to  pole 

75  to  KP*74*3 

U(LON)  =  U  wind  of  each  hemisphere  for  each  latitude 

at  each  sigma  layer 

V(LON)  =  V  wind  of  each  hemisphere  for  each  latitude 
at  each  sigma  layer 

Z(LON)  =  Z  height  of  each  hemisphere  for  each 
latitude  at  each  sigma  layer 

These  three  records  are  repeated  by  northern 
hemisphere  to  southern  hemisphere  for  37  latitudes 
(from  equator  to  pole)  for  KP=12  sigma  layers 

2.  FORTRAN  UNIT  2  -  Analyzed  residuals  at  regular  global  grids  (2.5°  X 

2.5°) 

(See  output  file  FORTRAN  UNIT  30  of  program 
C7DIXANL.) 
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OUTPUT  FILES: 


1.  FORTRAN  UNIT  10  -  The  reconstructed  global  grid  of  variables  U,  V  and 

temperature  at  updated  sigma  layers 


Record  #  Description 

1  to  KP 

KP=total  number 

of  sigma  layers  XT ( NLON , NLAT )  =  Temperature  on  regular  global  grid 

(2.5°  X  2.5°)  for  KP=12  sigma  layers, 
where 

NLON  =  144 
NLAT  =  73 


KP+1  to  2*KP  Ul(NLON.NLAT)  =  U  on  regular  global  grids 

(2.5°  X  2.5°)  for  KP=12  sigma 
layers 

2*KP+1  to  3*KP  VI (NLON, NLAT)  =  V  on  regular  global  grid 

(2.5°  X  2.5®)  for  KP=12  sigma 
layers 

3*KP+1  XQ( NLON, NLAT)  =  Surface  pressure  on  regular  global 

grid  (2.5®  X  2.5®) 


DESCRIPTION  OF  MAIN  SUBROUTINE: 


Designators 

CA 

CB 

CC 
CD 
CE 

Subroutine  FMTTRA  is  called  to  rearrange  the  data 
structure. 

Subroutine  PSZTOT  is  called  to  compute  update 
surface  pressure  and  to  convert  height  field  to 
temperature  field. 

Subroutine  TOSIG  is  called  to  update  variables  to 
new  sigma  layers  with  new  surface  pressure. 


Text 

Set  up  constants,  resolution,  sigma  structure  and 
grid  locations. 

Read  in  FGGE  3A  terrain  height  field  ¥ ( NLON , NLAT ) 
from  Tape  4. 

Get  first  guess  field  of  post-processing  from  Tape  1. 

Get  analyzed  field  from  Tape  2. 

Subroutine  SETSIG  is  called  to  set  up  sigma 
structure. 


PROGRAM: 


C9SBSPDATA 


ABSTRACT 

MAIM  SUBROUTINE: 
SUBROUTINE  CALLED : 


INPUT  FILES: 

1.  FORTRAN  UNIT  1 

2.  FORTRAN  UNIT  3 

Records  # 

1 

3.  FORTRAN  UNIT  A 

Records  # 

1  to  KP 

OUTPUT  FILES: 

1.  FORTRAN  UNIT  2 


This  program  converts  regular  global  grid  data  (2.5° 
X  2.5°)  to  spectral  data.  The  regular  grids  are 
linearly  interpolated  to  Gaussian  latitudes  first, 
then  the  field  is  transformed  to  the  spectral  domain. 

SPDATA 

LINTERP 

LEGSUM 

GAL 

LMN 

FFT1 

PMNS 

PRINT 

GAUSLAT 

POLY 


Regular  global  grid  (see  output  file  FORTRAN  UNIT  10 
of  program  C8GRIDF) 

Spectral  terrain  height  from  FGGE  3A  (rhomboidal 
truncation) 

Description 

ZMNl (MS , NS , 1)  =  Terrain  height  in  spectral  form, 

where  MS=31  and  NS=31  for  rhomboidal 
truncation  at  wavenumber  30 

Specific  humidity  in  spectral  form  from  FGGE  3A 
(rhomboidal  truncation) 

Description 

ZMNl (MS, NS, KP)  =  Specific  humidity  in  spectral  form 
for  rhomboidal  truncation  at 
wavenumber  30  with  KP=12  sigma 
layers 


Spectral  data  of  variables  vorticity,  divergence, 
temperature,  specific  humidity,  surface  pressure  and 
surface  geopotential;  in  coded  form  with  standard 
structure  (See  input  file  FORTRAN  UNIT  1  of  program 
C2P0STPTZ . ) 
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DBSCRIPTIOjj  OF  MAIM 
Designators 
CA 

CB 

CC 

CD 

CB 


PROGRAM: 


Text 

Set  up  constants,  resolution,  Gaussian  quadrature 
and  Fourier  transform. 

U  and  V  grids  are  linearly  interpolated  to  Gaussian 
latitudes  and  are  converted  to  vorticity  and 
divergence  in  spectral  form. 

Compute  spectral  components  for  temperature . 

Merge  specific  humidity  data. 


CF 


Compute  spectral  components  for  surface  pressure. 
Merge  surface  geopotential  data. 


/  / 


PROGRAM: 


ABSTRACT: 


MAIM  SUBROUTINE: 


This  program  is  modified  from  NMC  global  forecast 
model  to  perform  global  normal  mode  initialisation. 

GMAIN  -  This  main  subroutine  is  set  up  to  get  initial 
parameters,  coefficients  constants,  model  structures 
and  file  assignments  for  model  processing.  (See  GSM 
program  documentation  for  details^*).) 

GNMINI  -  Subroutine  is  called  by  subroutine  GMA1M 
to  apply  Machenhauer  method  of  nonlinear  normal  mode 
initialization  for  global  data  set.  The  normal 
modes  with  vertical  mode  less  than  MODS  and  period 
less  than  PERCUT  are  adjusted  to  produce 
approximately  zero  tendency. 


SUBROUTINES  CALLED: 


1.  GMAIN: 


SETSIG 

AMHMTM 

GLATS 

EPSLON 

PRMFLD 

GRDLNF 

GSSTCD 

GNMINI 

GWRITE 

RMS 

BMCM 

GLOOP 

TEND 

VERTIC 

PRIMES 

H0RIZ1 


The  details  of  all  these  subroutines  are  described  in  GSM  program 

_ .  (A)  _  Lt. _  _  C _ 1  2 _ _ 1 -J  _  i.lLt.ll L  i I 


2.  GNMINI: 


documentation . 


The  theory  of  nonlinear  normal  mode  initialization  is 


discussed  in  Ballish,  1980v  and  generalization  of  the  scheme  is 
described  in  Gerlach,  1983. ^ 


4.  GSM  Program  Documentation,  NMC/AWS,  1983. 

5.  Gerlach,  A.  M.  (ed.),  1983:  Objective  Analysis  and  Prediction 

Techniques  -  1983.  APGL-TR-83-0333,  Contract  F19628-82-C-0023,  Systems  and 

Applied  Sciences  Corporation,  ADA142441. 
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INPUT  FILES: 


1.  FORTRAN  UNIT  5  -  A  record  containing  28  integers  in  the  format  2815. 

These  integers  fill  the  array  HUM 

The  details  of  the  array  HUM  are  described  in  GSM  program 

(4) 

documentation.  The  inputs  NUM(18)  and  NUM(19)  are  used  in  the 
initialization  scheme. 

NUMC18)  =  MODS;  Number  of  vertical  modes  to  initialize 

NUM(19)  =  NITER:  Number  of  iterations  for  Machenhauer  normal  mode 
initialization 


2.  FORTRAN  UNIT  16  -  Inputs  for  sea  surface  temperature  and  drag 

coefficients 


Record 


Description 


Sea  surface  temperature  with  dimension  (#LON,#LATG) 
ordered  from  equator  to  poles;  sea  surface 
temperature  less  than  269.95  is  considered  over  land 

Drag  coefficients  with  dimension  (#LON,#LATG) 
ordered  from  equator  to  poles, 

where  #LON  =  Number  c c  longitudes  for  each  latitude 
#LATG  =  Number  of  Gaussian  latitudes 


Both  records  are  reordered  in  subroutine  ASSTCD  from  north  pole  to 
south  pole. 

3.  FORTRAN  UNIT  18  -  Spectral  data  input  (unformatted)  (See  FORTRAN  UNIT 

1  in  program  C11P0STPTZ.) 

4.  FORTRAN  UNIT  80  -  Inputs  for  eigenvalues  and  eigenvectors  from  the 

normal  mode  computation 


Record  # 


Description 


EIGG , EIGGT , GH , D0TPR0 , TO 

EIGG(#KP,#KP)  =  The  sorted  unit  eigenvectors  of 
matrix  G 

EIGGT (#KP,#KP)  =  The  sorted  unit  eigenvectors  of 
matrix  G^ 

GM(fkp)  =  The  negative  of  the  eigenvalues  of 
matrix  G 

D0TPR0 (#KP)  =  The  reciprocal  of  the  inner  product 
of  the  unit  eigenvectors  of  G  and  G^ 


TO(#KP)  s  A  vector  of  basic  state  temperatures  used 
in  the  generation  of  the  normal  modes, 

where  #KP  =  Number  of  vertical  layers  or  modes 

2  to  #KP*#JCAP1+1  PER.G 

where  #KP  =  number  PER(JG)  =  Periods  of  gravity  modes  for  each  zonal 
of  vertical  modes  wave  and  each  vertical  mode, 

#JCAPl=Zonal  wave  where  JG  =  Number  of  gravity  modes 

number  of  rhom- 

boidal  truncations  G(JG,NAS)  =  Corresponding  eigenvectors, 

fJCAP  +  1. 

where  NAS  =  Vector  size  of  symmetric  or 
antisymmetric 

This  pattern  is  repeated  for  fJCAPl  zonal  wavenumber 
and  #KP  vertical  modes. 

OUTPUT  PILES: 

1.  FORTRAN  UNIT  19  -  Normal  mode  initialized  spectral  coefficients  for 

all  model  variables  (unformatted)  (See  FORTRAN 
UNIT  1  in  program  C11POSTPTZ.) 

DESCRIPTION  OF  MAIN  SUBROUTINES: 

See  GSM  program  documentation^* )  for  details. 


PROGRAM:  C12RPOSTP 

ABSTRACT;  This  program  is  set  up  to  compute  the  update 

residuals  at  observation  sites.  The  normal  mode 
initialized  field  is  post-processed  at  observation 
sites  and  combined  with  original  observation  data  to 
construct  update  residuals.  Sela’s  method^)  is 
used  to  convert  temperature  residuals  to  height 
residuals . 

MAIM  SUBROUT  liTK:  RPOSTP 

SUBROUT IWKS  CALLED;  LEGSUM 

FFTGP 

PMMS 

BSCST 

SPTOCP 

UMVM 

SIGMMC 

GETDZ 

DTTODZ 


INPUT  FILES; 

1.  FORTRAN  UNIT  1  -  Normal  mode  initialized  spectral  data  (unformatted) 

(See  FORTRAN  UNIT  1  in  program  C11POSTPTZ.) 

2.  FORTRAN  UNIT  10  -  Observations  of  U,  V  wind  and  temperature  with  equal 

number  of  points  on  each  layer  (See  FORTRAN  UNIT  15 
in  program  C4RFMPR0  for  details.) 

3.  FORTRAN  UNIT  20  -  A  matrix  and  C  matrix  for  Sela's  hydrostatic 

equation  (See  FORTRAN  UNIT  10  in  program  C4RMFPR0 . ) 

OUTPUT  FILES: 

1.  FORTRAN  UNIT  7  -  Updated  residuals  at  observation  sites 
Record  #  Description 

1  N  =  Number  of  observations  at  layer  1  for  U 

2  X(I)  =  Latitudes  for  U 

3  X(I)  =  Longitudes  for  U 

4  X(I)  =  U  residuals  at  observation  sites  on  sigma 

layer  1 

5  N=  Number  of  observations  on  layer  1  for  V 


6 


X(I)  =  Latitudes  for  V 


7 


X(I)  =  Longitudes  for  V 


8  X(I)  =  V  residuals  at  observation  sites  on 

sigma  layer  1 

9  H  =  Number  of  observations  at  layer  1  for  composite 

variable  P 

10  X(I)  =  Latitudes  for  P 

11  X(I)  =  Longitudes  for  P 

12  x(l)  =  Residuals  of  composite  variable  P  at 

observation  sites  on  sigma  layer  1 

These  twelve  records  are  repeated  for  all  sigma  layers. 

DESCRIPTION  OF  MAIM  SUBROUTINE: 


Desiznators 

Text 

CA 

Set  up  basic  constants  and  resolutions. 

CB 

Read  in  initialized  data. 

CC 

Input  for  observation  locations. 

CD 

Post- process  initialized  data  at 
and  store  first  guess  on  scratch 

observation  sites 
files. 

CE 

Check  date  and  compute  residuals 
T. 

of  variables  U,  V, 

Call  subroutine  GETDZ  to  get  residuals  of  height 
then  convert  to  residuals  of  composite  variable. 
The  results  are  rearranged  and  written  to  output 
file. 


CF 


PROGRAM: 

ABSTRACT: 


CRUMP 

This  program  is  designed  to  calculate  the  normal 
mode  functions  for  the  AFGL  model  to  be  used  as 
basic  functions  in  the  analysis  of  velocities, 
heights,  and  surface  pressure.  The  coefficients  of 
vorticity,  divergence  and  the  composite  variable  P 
are  read  in  from  PORTRAIT  UNIT  1  for  the  S=0  case  and 
from  PORTRAIT  UNIT  2  for  S>0  case.  The  normal  mode 
functions  are  stored  in  PORTRAIT  UNIT  5.  Both  Rossby 
wave  and  gravitational  modes  with  periods  greater 
than  48  hours  are  included. 


CRNMP 


SUBROUTINES  CALLED: 


MODSUM 

PLOT 

PNMS 

WROUT 


INPUT  PILES: 

1.  FORTRAN  UNIT  1  -  The  coefficients  of  vorticity,  divergence  and 

composite  variable  P  for  S=0 


Record 


Description 


GH(IM)  »  Height  of  vertical  mode 


2  to  3 


(PERU ,  IAS) ,  I-NP1  ,NPST)  , 

( (G(II, JJ.IAS) , II=NFI,NFST) , 
JJal.NMER) 


PBR(I.IAS)  a  Periods  for  Rossby  waves  only, 
where  IAS  is  index  of  symmetric  and  antisymmetric 
modes 

IASal  for  symmetric  and 
I AS =2  for  antisymmetric 

G(II , JJ , IAS)  a  Coefficients  for  Rossby  waves  only, 
where  HMER  is  the  number  of  coefficients  found  for 
each  eigenvector  Which  includes  total  of  divergence, 
vorticity  and  composite  variable  P 

The  records  1  to  3  are  repeated  over  all  the  desired  vertical  modes 
IMODBS. 


2.  FORTRAN  UNIT  2  -  The  coefficients  of  vorticity,  divergence  and 

composite  variable  P  for  S>0;  both  Rossby  and 
gravitational  modes  (whose  period  exceeds  PERCUT) 
are  included 
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1 


( (BLANK (I , J) , I=ILEV) , J=1 , I LEV) , 
((BLAMKU.J)  ,1-l.ILBV) , J-l.ILEV) , 
(GH(I),I=1,ILEV) 


GH(I)  =  Height  of  vertical  modes 

2  to  3  (PER(I , IAS) ,1=1 ,NF) , ( (G(II, JJ.IAS) , 

11=1 ,NF) , JJ=NMER) 

PER(I.IAS)  =  Periods  for  both  Rossby  waves  and 
gravity  waves 

G(II,JJ,IAS)  =  Coefficients  for  both  Rossby  waves 
and  gravity  waves 

The  record  2  to  3  are  repeated  over  all  the  zonal  waves  and  all  the 
vertical  modes. 

OUTPUT  FILES: 

1.  FORTRAN  UNIT  5  -  The  normal  mode  functions  for  every  1  degree 

latitude  interval 

(See  input  file  FORTRAN  UNIT  2  in  program  C6DIZEXPN0 
for  details.) 

2.  FORTRAN  UNIT  3  -  Total  number  of  gravity  waves  and  Rossby  waves  to  be 

stored  for  each  zonal  wave  at  each  vertical  mode 


Record 


Description 


1  to  IMODES*NSWV  LALL  =  Number  of  frequency  indices  for  each  zonal 

where  IMODES  is  wave  and  each  vertical  mode 

the  total  number  of 

vertical  modes,  NSWV 

is  total  number  of 

zonal  waves 


DESCRIPTION  OF  MAIN  SUBROUTINE: 

Designators  Text 

CA  Read  in  all  vertical  modes. 

CB  Write  geopotential  modal  heights  to  output  file. 

CC  Read  in  periods  and  coefficients. 

CD  Compute  Legendre  functions  for  all  latitudes  and 

specific  S. 


Das  Unitors  Text 

gg  Compute  normal  mode  functions  for  Rossby  modes  and 

call  subroutine  WROUT  to  write  to  output  file. 


CF 


Compute  normal  mode  functions  for  gravitational 
modes  and  write  to  output  file. 


PROGRAM: 


ZEROS 


ABSTRACT :  This  program  is  set  up  to  calculate  coefficients  of 

vorticity,  divergence,  and  composite  variable  for 
zonal  wave  S*0. 

MAIM  SUBROUTINE:  ZEROS 

SUBROUTINE  CALLED:  Mone 

OUTPUT  FILES: 

1.  FORTRAN  UNIT  1  -  The  coefficients  of  vorticity,  divergence,  and 

composite  variable  for  S=0 

(See  input  file  FORTRAN  UNIT  1  in  program 
NMFUNCTIONS . ) 


DESCRIPTION  OF  MAIM  SUBROUTIMB: 

Designators  Text 

CA  Input  of  vertical  modes  from  data  statement  and 

output  of  these  to  FORTRAN  UNIT  1. 

CB  Loop  over  all  the  Rossby  frequencies  and  fill  in 

sparse  vector.  Set  divergence  coefficient  to  zero. 

CC  Perform  Gramm-Schmidt  procedure. 

CD  Compute  coefficients  of  vorticity  and  composite 

variable,  then  write  to  output  file. 


*  -  •  '  »  •  < 
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PROGRAM: 


ABSTRACT: 


BIGCOM 

RBIGCOM 

These  programs  are  set  up  to  compute  the  normal  mode 
coefficients  for  use  in  nonlinear  normal  mode 
initialization.  The  output  of  program  EIGCOM 
contains  only  gravity  modes.  Both  Rossby  modes  and 
gravity  modes  are  included  in  REIGCOM  output. 


The  details  of  these  programs  are  described  in  GSM  program 
documentation. The  basic  theory  is  in  Ballish  (1980). 


IV.  DESCRIPTION  OF  SUBPROGRAMS 

The  details  for  those  subroutines  called  by  the  main  programs  are 
discussed  in  this  section.  It  is  written  in  the  following  form  for  each 
subroutine: 

SUBROUTINE :  (Name  of  subprogram) 

ABSTRACT :  (Abstract  of  the  subroutine) 

SUBROUTINES  CALLED:  (Name  of  subroutines  called  by  this  subroutine) 
INPUT  VARIABLES:  (Variables  used  as  input  arguments) 

OUTPUT  VARIABLES:  (Variables  used  as  output  arguments) 

ARGUMENT  LIST:  (Variables  in  argument  list) 


SUBROUTINE : 
ABSTRACT: 

SUBROUTINES  CALLED: 
INPUT  VARIABLES: 
Symbol 
DTCK) 

KP 

CH 

N 

OUTPUT  VARIABLES: 
Symbol 
DX(k) 


DTTODZ 

This  subroutine  is  set  up  to  get  height  residuals 
from  temperature  residuals  by  using  Sela's  method. 

None 

Meaning 

Temperature  residual  profile  with  dimension  KP 
Model  vertical  resolution 
C  Matrix  from  Sela's  hydrostatic  equation 
Number  of  observations 

Meaning 

Variable  to  store  height  residuals  with  k-1  at 
bottom 


DZ(k) 


Computed  height  residuals  from  Sela's  hydrostatic 
relation,  where  k=l  at  top 


SUBROUTIWE: 


FLATZT 


ABSTRACT:  This  subroutine  is  set  up  to  apply  Flattery’s 

algorithm^1*  to  get  layer  temperatures  from  layer 
heights.  The  method  involves  the  construction  of 
thickness  temperatures  from  heights.  Using  these 
thickness  temperatures,  a  set  of  equations  is 
obtained  by  least-squares  approximation  that  yield 
temperatures  which  minimize  the  errors  of  the 
constraining  equations. 

SUBROUTINE  CALLED:  IMINV 


INPUT  VARIABLES: 
Symbol 
ZL(MLPl) 
SL(HLPl) 
MLP1 
IWRT 

OUTPUT  VARIABLES: 
Symbol 
TL(MLPL) 


Meaning 
Layer  heights 
Model  sigma  layers 
Number  of  layers 
FLAG  for  printer 

Meaning 

Layer  temperatures 


Description: 

DO  Loop  20  constructs  the  reciprocal  of  the  difference  across  the 
layers  of  the  logarithm  of  pressure. 

DO  Loop  30  sets  up  coefficients  in  one  set  of  constraining  equations 

IX)  Loop  40  initializes  to  zero  all  the  elements  of  the  matrix  A 
transpose,  AT(AT) . 

DO  Loop  50  and  60  set  up  the  coefficients  of  matrix  AT. 

DO  Loop  80  initializes  to  zero  the  elements  of  the  matrix  product 
ATA  or  ATA,  and  then  calculates  the  matrix  product. 


The  thickness  temperatures  TBAR  are  computed  from  the  input  heights 
and  then  the  right  hand  sides  of  the  constraining  equations  are 
computed  in  UNK. 


The  right  hand  side  of  the  equations  AT  *  U  is  multiplied  by  AT. 

The  solution  is  obtained  by  multiplying  the  inverse  of  ATA  by 
A?U  and  scaled  by  multiplying  by  SCLT. 


subroutine 


FHTTRA 


ABSTRACT; 

SUBROUTINE  CALLED 
INPUT  VARIABLES: 
Symbol 
NFl 
NF2 
A 

NIN 

NOUT 

NLON 

OUTPUT  VARIABLES: 
Symbol 
B 


This  subroutine  is  called  to  rearrange  the  data  by 
changing  dimension. 

None 

Meaning 

Original  unit  number 
Rewrite  unit  number 
Original  array 
Original  dimension 
Redefine  dimension 
Fixed  dimension 

Meaning 

Rearranged  array 


FUVCH 


ABSTRACT: 

SUBROUTINE  CALLED 
INPUT  VARIABLE: 
Symbol 
PHI 
IK 
XI 
X2 
I  FUN 

IAS 

OUTPUT  VARIABLE: 
Symbol 
XRET 


This  subroutine  is  set  up  to  obtain  basic  function 
at  given  latitude  frost  linear  interpolation  of  basic 
function  table.  The  table  is  constructed  with  every 
one  degree  latitude  intervals. 

None 


Meaning 

Given  latitude 

Index  for  interpolation 

Lower  bound  basic  function 

Upper  bound  basic  function 

Function  indicator 

IFUM=1  for  U 
IFUM=2  for  V 
IFUN=3  for  Z 

Index  for  symmetric  or  antisymmetric  modes 

IAS=1  for  symmetric 
IAS=2  for  antisymmetric 

Meaning 

Basic  function  at  given  latitude 


SUBROUT  IITK: 


GETDZ 


ABSTRACT; 

SUBROUTINE  CALLED: 
ARGUMENT  LIST: 
Symbol 
X 

RES 

NALL 

CM 

DX 

DZ 

QX 


This  subroutine  is  called  to  compute  residuals  for 
composite  variable  and  construct  the  output  file. 

DTTODZ 

Meaning 

Variable  array  for  input  or  output 
Residual  array  of  composite  variable 
Number  of  data  points  for  each  sigma  layer 
C  Matrix  from  Sela's  hydrostatic  equation 
Height  residuals  profile  from  bottom  up 
Height  residuals  profile  from  top  down 
Array  of  surface  pressure  residuals 


I 


SCBgOUTIHg: 


GBTOBS 


ABSTRACT: 

SUBROUTIMBS  CAL 
ARCUMRMT  LIST: 
Symbol 
OBS 

KTOTAL 

ALAT 
A  LON 
CCS 
NFUN 
KFUN 

(CALL 


This  subroutine  is  called  to  read  in  normal  mode 
data  projections  from  input  file. 

Bone 


Meaning 

Array  for  variable  input 

Humber  of  data  points  for  each  variable  at  each 
vertical  mode 

Array  of  data  latitudes 

Array  of  data  longitudes 

Array  of  exponentials  at  data  longitudes 

Number  of  variables 

Number  of  data  points  for  each  variable  at 
each  vertical  mode 

Total  number  of  data  points  for  all  the  variables 
at  each  vertical  mode 


SUBROUT I MB: 


GBTPS 


ABSTRACT: 

SUBROUTINES  CALLED: 

IMPUT  VARIABLE: 
Symbol 
PHI (LAYER) 
LAYER 
KP 
PO 
SL 
P 
PP 
N 


This  subroutine  is  set  up  to  obtain  updated  surface 
pressures  by  using  quadratic  relationship. 

Mone 

Meaning 

Height  profile  with  dimension  LAYER 

Number  of  layers  included  surface 

Number  of  sigma  layers 

Model  surface  pressure 

Sigma  value  at  layers 

Array  of  log  pressure 

Array  of  pressure 

Number  of  data  points 


OUTPUT  VARIABLE: 

Symbol  Meaning 

DQ  Surface  pressure  residual 


PN 


Updated  surface  pressure 


SUBROUTINE: 


GBTTH 


ABSTRACT :  This  subroutine  is  set  up  to  obtain  the  appropriate 

horizontal  normal  mode  functions  for  analysis. 

SUBROUT IRES  CALLED:  PUNCH 

INPUT  VARIABLES: 


Symbol 

Mean ins 

PH 

Normal  mode  functions  at 
intervals 

every  one  degree  latitude 

RANLAT 

Array  of  data  latitudes 

KTOTAL 

Number  of  data  points  at 
all  the  variables 

each  sigma  layer  for 

KFUN 

Number  of  data  points  at 
variable 

each  sigma  layer  for  each 

NLAT 

Not  in  use 

NPUN 

Number  of  variables 

IAS 

Index  for  symmetric  or  antisymmetric 

OUTPUT  VARIABLE: 

Symbol 

Mean inn 

THET 


Array  of  normal  mode  functions  at  observation  sites 


SUBROUT  IlfE : 


GBTZBRO 


ABSTRACT: 

SUBROUTIireS  CALLED: 
ARGUMKMT  LIST: 
Symbol 
X 

OBSLAT 

OBSLON 

NALL 

KP 

Nl 

N2 

IRES 


QX 


This  subroutine  inserts  zero  residuals  at  the 
corners  of  a  grid  box  if  no  observation  is  found  in 
that  box. 

None 

Mean ins 
Array  of  data 
Array  of  latitude 
Array  of  longitude 

Number  of  data  points  at  each  sigma  layer 

Number  of  sigma  layers 

Original  unit  number  of  data  file 

Final  unit  number  of  data  file 

FLAG  for  residuals 

IRES-0  is  for  observation 
IRES=1  is  for  residual 

Array  of  updated  surface  pressure 
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SUBROUT  IlfK: 


IMIHV 


ABSTRACT: 


This  subroutine  is  called  to  compute  inverse  of 
matrix.  The  standard  Gauss-Jordan  method  is  used. 


SUBROUTIHBS  CALLED: 
IMPUT  VARIABLES: 
Symbol 
A(»,H) 

N 

L(H) 

M(N) 

OUTPUT  VARIABLES: 
Symbol 
A(H.N) 


Hone 

Mean ins 

Input  matrix  of  order  M  to  be  inverted;  destroyed 
computation  and  replaced  by  resultant  inverse 

Order  of  A 

Work  vector  of  length  II 
Work  vector  of  length  V 

Meaning 

Inverse  of  input  matrix  A 

Determinant  of  A  inverse 
If  D=0.0,  A  is  singular 


D 


SUBROUT I MB: 


PSZTOT 


ABSTRACT: 


SUBROUTINES  CALLED: 
IMPUT  VARIABLES: 
Symbol 
PHI (LAYER) 
LAYER 
KS 
PO 
SL 
PL 
P 

PP 

IWRT 

OUTPUT  VARIABLES: 
Symbol 
T(KS) 


This  subroutine  is  set  up  to  compute  updated  surface 
pressure  by  using  quadratic  relationship. d)  The 
updated  surface  pressure  is  the  input  for  routine 
FLATZT  to  convert  height  profiles  to  temperature 
profiles. 

FLATZT 


Meaning 

Height  profile 

Number  of  layers  including  surface  layer 

Number  of  sigma  layers 

New  surface  pressure 

Sigma  values  at  layers 

Update  pressures  at  layers 

Array  of  log  pressure 

Array  of  pressure 

FLAG  for  printer 

Meaning 

Updated  temperature  profile 
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SUBROUTINE: 

ABSTRACT: 


SETSIG 


SUBROUTINES  CALLED: 
Symbol 
KP 
KPP1 
KPN1 
MK 
MKM1 

OUTPUT  VARIABLES: 
Symbol 
DEL(KP) 
SI(KPPl) 

SL(KP) 

CI(KPPl) 

CL(KP) 


This  subroutine  is  called  to  set  up  vertical  sigma 
structure. 

None 

Meaning 

Number  of  sigma  layers 
KPPl  =  KP  +  1 
KPM1  =  KP  -  1 
Not  in  use 
Not  in  use 

Meaning 

Sigma  spacing  for  layers 
1.0-CI  at  levels 
Sigma  values  at  layers 
Sigma  values  at  levels 
1.0- SL  at  layers 
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SIGNMC 


SUBROUTINE: 

ABSTRACT : 

SUBROUTINES  CALLED: 
INPUT  VARIABLES: 
Symbol 
U 
V 
T 

GZ 

IHEM 

OUTPUT  VARIABLE: 
Symbol 
US 
VS 
TS 
GZS 
DELNMC 
SINMC 


This  subroutine  is  called  to  change  the  variables 
from  AFGL  sigma  structure  to  NMC  sigma  structure 
(index  increasing  upward). 

None 

Meaning 

Array  of  U 

Array  of  V 

Array  of  temperature 

Surface  geopotential 

IHEM=1  for  northern  hemisphere 
IHEM=2  for  southern  hemisphere 

Meaning 

Array  of  U  from  surface  upward  (NMC) 

Array  of  V  from  surface  upward  (NMC) 

Array  of  temperature  from  surface  upward  (NMC) 

Surface  height  in  meters 

DEL  from  surface  upward  (NMC) 

SI  from  surface  upward  (NMC) 

SL  from  surface  upward  (NMC) 


SLNMC 


SUBROUTINE: 


TOSIG 


ABSTRACT: 

SUBROUTINES  CALLED: 
ARGUMENT  LIST: 
Symbol 
XP 
XS 
PS 

NLON 

KP 

MK 

MKM1 

VLNP 

DLNP 

SL 

PP 

PL 

ALNP 

PM 


This  subroutine  is  called  to  interpolate  the 
variables  to  new  layers  with  updated  surface 
pressure. 

Hone 

Meaning 

Array  of  variable  at  original  vertical  layers 
Array  of  variable  at  updated  vertical  layers 
Updated  surface  pressure 
Longitudes 

Number  of  new  layers 
Number  of  old  layers 
MKMl  s  MK  -  1 

Mean  of  layer  log  pressures 
Difference  of  layer  log  pressure 
Sigma  value  of  vertical  structure 
Original  layer  pressures 
Updated  layer  pressures 
log  layer  pressures 
original  surface  pressure 


SUBROUTINE : 


TTOZ 


P 


i 

i 


ABSTRACT: 


SUBROUTINES  CALLED: 
INPUT  VARIABLES: 
Symbol 
TT 
ZZ 
DZ 
KP 

KPP1 

CM 

AM 

S 

ZSTAR 

N 


This  subroutine  is  set  up  to  compute  layer  heights 
from  layer  temperatures  by  using  Sela's  hydrostatic 
relations. ^ 

None 

Meaning 

Layer  temperatures  from  surface  upward 
Layer  temperatures  from  top  down 
Layer  heights  from  surface  upward 
Number  of  layers 
KPP1  =  KP  +  1 

C  Matrix  of  Sela's  hydrostatic  equation 
A  matrix  of  Sela's  hydrostatic  equation 
Vector  for  surface  height 
Surface  height 
Nth  data  point 


k* 
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.V.wV  .'W  • 

mUmUiAA 


\v.yv  ■ 


SUBROUTINE: 


VERT 


ABSTRACT: 

SUBROUTINES  CALLED: 
ARGUMENT  LIST: 
Symbol 
OBS 
MOBS 

KOBS 
IMODES 
I  FUN 


This  subroutine  is  set  up  to  compute  the  vertical 
projection  of  the  variable  found  in  the  array  OBS 
and  to  store  it  in  OBS  with  layers  changed  to  modes. 

None 


Meaning 

Array  of  variable 

Maximum  dimension  for  array  OBS  to  be  used  in  this 
subroutine 

Number  of  data  points  for  the  variable 

Number  of  vertical  modes  used 

Variable  indicator 

IFUN=1  for  U 
IFUN=2  for  V 
IFUN=3  for  Z  (height) 


NLEV 


Number  of  layers 


i>4.--»«iy^iwi  j 


'  f*.  4 


SUBROUTINES : 

AFGL  IN-HOUSE  SUBROUTINES: 

BSCST 

FFT1 

GAUSLAT 

GQL 

LEGSUM 

LINTERP 

LMN 

PMNS 

POLY 

PRIN 

SPTOGP 

UMVM 

01  PACKAGE  SUBROUTINES:  (PROGRAM  C3ASAPRES ) 


ASAPl 

MAST0R1 

MAST0R4 

QCALC 

OBTMZL 

MAST0R2 

MAST0R6 

CALCRES 

FLAGS 

FG 

SPTOGPW 

LEGSUM 

PMNS 

UMVM 

CQCV 

PTOSIG 

TLAZLE 

FFTGP 

SPTOGPX 


Not  all  subroutines  from  AFGL  in-house  package  or  01  package  are 
described  in  this  document. 


V.  SUMMARY  OF  RUN  STREAM 


The  run  stream  for  the  experiment  is  summarized  in  Table  1.  This 
table  is  an  outline  from  which  a  procedure  file  may  be  developed  to 
implement  the  run  stream.  Corresponding  COSMOS  procedure  files  for  the 
AFWL  Cray-1  computer  system  are  displayed  in  Figure  1  and  Figure  2. 
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•  *  v"  v'  «  '  •  .V  . 


•  select  prinUo*fctlo9,sovef>holtenj,tesk«4« 

2  •interrupt  on  softuoreerr  to  stp 

3  nfilo  s»ne«cositr,end«* 

J  **elect  pn#t lo*  1  jfc i<4> (tivtf*roiu It*, titM tfl 
5  •interrupt  on  ooftuoroorro  to  M 

*  *if  <4>  .eg.  1  then  fo  to  ski 
2  *ipU  /  2 

8  • destroy  tepeS  f*iat 
8  •twitch  tape2  fjint 

10  «n»es  store  fjint t/results/fjint<4> 

11  osuitch  terSI  topeS 

12  *iet  two  ■“20* 

13  oif  <1>  .90.  It  then  lot  too  ■  “2* 

14  *«im  jot  tope2sf99o\too\<1><2> 

15  ofile  none* input 

Id  fdtn  idnte>< l > , it ise><2>  fond 
17  txitp  /  It 

16  *ntt  store  tepe64:/resolts/nsop<4> 

19  *switch  topeS  terSI 

20  *90  to  sk2 

21  oskls  noss  jet  initi/results/init1?l2z 

22  *90  to  13 

23  *sk2i  destroy  topel  tapelt  tope 15  input 

24  osuitch  toped*  topel 

25  *switch  coot  topelO 

26  oxc4r  /  2 

2/  *noss  store  output t/rosults/c4op<4> 

28  enoss  store  topoSt/rosolts/zeroo«<4>  topel Si /results /otseq<4> 

29  *switch  output  c4op<4> 

30  ojive  c4op<4>  I486  ood 

31  edestroy  topel 
12  elet  itr*3 

S3  •snitch  topelO  coot 

34  «suitch  tope t 5  otsoq 

35  •switch  tope5  topel 
34  «suitch  reij  topeS 
37  oxver 

31  *11000  store  output i/results/verop\itr\<4> 

19  osuitch  output  verop\itr\<4> 

40  ojivo  verep\itr\<4>  1414  ood 

41  odestroy  topel  topo2  topel  topo20 

42  osuitch  topelO  topel 

43  osuitch  uup  topo2 

44  osuitch  loll  topeS 

45  ondin  /  19 

44  ones*  store  output i/results/diuop\ltrV<4> 

47  osuitch  output  divop\ltr\<4> 

48  ojivo  diuop\itr\<4>  I486  end 

49  enoss  store  tope20(/rosslts/rjwi\itr\<1><2><3> 

50  ocopy  tope24  tp28\itr\<4> 

Figure  1.  COSMOS  procedure  file  for  the  experiment  run 
without  iteration 


I 

1 

St  *9iv«  tp20\itr\<4>  1144  *M 

52  M#«trey  Up* 4 

53  MMitch  tap#2  uvp 

54  ••vitch  MV 17  tap*  4 

I  55  Mc7d  /  2 

54  •••••  •tor*  output i /r*»Ml t»/c 7op\itr\<4> 

57  •«••«  »tort  t*pt30t/rMult«/rfr\itr\<1><2Kl> 

SI  ••uitch  output  c?op\itr\<4> 

59  •fiu*  c7op\itr\<4>  1484  on! 

40  *ftvitch  top*5  rtij 

41  ••uitch  top*3  1*11 

42  *«Mitch  tapo4  uv37 

43  •4*«troy  tapol  tip#?  tip«10 

44  Muitch  fjiat  tapol 

45  ••Mitch  tap*30  tip* 2 
44  MMitch  torht  topo4 
47  ouclf  /  2 

41  *•!••  «tor*  output i/r*OMlt«Sc8op\itr\<4> 

49  •*•••  «tor*  tap*  1 0 1  /r*«M  1  to/fr id\ t tr\<1  K2XI> 

70  Moitch  output  c8op\itr\<4> 

71  «9iv*  cBop\itr\<4>  1484  ond 

72  «d*«troy  tapat  tap*2  tap*3 

73  Muitch  tap*4  t*rht 

74  «suitch  t*r31  tap*3 

73  ••uitch  sp<4>  tap*4 

74  ••uitch  tap*10  tapol 

77  *xc9»  /  4 

78  ***••  store  output t/r*«ult*/c9\itr\<4> 

79  ••Mitch  output  c9op\itr\<4> 

00  Miv*  c9*p\itr\<4>  1484  ood 
•1  «SMitch  tap*4  sp<4> 

•2  ••Mitch  tap* 3  tvrll 

13  •*»»  «tor*  tap#2i/ro«ult«/c9out\itr\<1X2X3> 

•4  «d*«troy  tap*5  tap* 14  tap*1l  tap*00  tapol 9 

IS  ••Mitch  mum  tapes 

04  ••Mitch  udp  tap# 14 

17  «»Mitch  tap*2  taptll 

IB  ••uitch  oif  tap*IO 

•9  •xpa*  /  4 

90  Milt  «t*r*  tap*19s/r*«MUs/iBlt\itr\<1X2><3> 

91  •«••«  •tor#  t*p*4i/MMltt79«out\itr\<IX2XI> 

92  ••uitch  tap*4  t4\ttr\<4> 

93  •9iv*  t4\itr\<4>  1484  t.  oad 

94  ••Mitch  t4\itr\<4>  tap* 4 

95  MMitch  tap#3  nun 
94  «tuitch  tap#80  *if 
97  MMitch  tap* 14  udg 

91  Han  «t*r#  tap*18s/r*ftults/«pd\itr\<1X2Xl> 

99  ••Mitch  tap#19  iait 

100  M##tr*y  tapcl  tap#2  tap* 10  output  tap* 20 

Figure  1.  COSMOS  procedure  file  for  the  experiment  run 
without  iteration  (continued) 
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10t  *13:  twitch  »d$2  tip*  2 

102  edtfttroy  ttptl 

103  ttwitch  ml  ttptl 

104  tdtttroy  tip* 2  ttptl 

105  *xfct  /  • 

104  tdtttroy  tip* 8 

107  *t«itch  top*  2  tdj2 

108  tltt  vd»<1> 

100  «ltt  vt  ■  -00" 

110  tif  <2>  .tq.  "00"  thtt  ltt  vt«"12" 

111  tif  vt  .tq.  “00“  ihtn  ltt  vd*vd*1 

112  *mm  tttrt  t*p*7i/rttultt//cUvd\\vt\ 

111  ttoitch  t»p*7  fct\vd\\vt\ 

114  tjivt  fct\vdN\vt\  1844  h.  ttd 

115  tdtttroy  ttptl 

111  tjo  tt  ltt 
117  *14t  ft  to 

111  tlttt  twitch  fct\vd\\vt\  ttptl 
Ilf  f 

120  tnttt 

121  fet  ff  itrll  tttrrtttllc  MMittptS  wdgttaptU  fnn  tifittptBON 

122  «d92ittpt2 

123  dtftult  dir"/htltotf 

124  §tt  coiticnitrixt  rtl|trtisl1i12  ovpitvpntfutl2  ltlltltllotfil2l 

125  uv37iuvpnnf371it32  ttrht :t*rht144  tp1ttp1712z  tp2ttp1800z  tp3ltp1812zl 
124  tp4itp1900z  tp5ttp1912z  tp4ttp21000z 

127  tot  dir*/txtcutt 

128  ftt  xdix  xfct  xc4r  xvtr  xc7d  xc8f  xc9t  xptt  xttp  xfim 

129  ttd 

1M  tltt  lot  •  7 

131  tltt  ft»0 

132  «lot  Jd  -  • 

133  tltt*  ltt  Jd  •  Jd«  1 

134  tif  ft  thtn  lot  it  •  *00" 

135  tif  .tot.  ft  thtt  lot  ft  ■  ”12" 

134  •cttntt  i*cttitr  tith  \tnc\,\*t\f"z*,\Jd\ 

137  tif  .tot.  ft  thtt  lot  itc  •  itc  ♦  1 
111  tltt  ft  •  .tot.  ft 

139  *if  Jd  .ft.  5  thtn  ft  to 

140  tft  to  lot 

141  o«tpi  ft  to 


Figure  1.  COSMOS  procedure  file  for  the  experiment  run 
without  iteration  (continued) 


1«*SEl|CT  PVlNTLQ*CYCllC,SVVEF*HAlTUNG,TASK*40 
2»*lWTCRfttfPT  ON  lerTMKEH  10  STP 
l**FliC  M*M£bC0S!1K,CMI>? 

4**S£L£C  T  PI IHTlO*tSI  Tft<4>  , SAVEF-AESULTS , TASK-10 
5»* INTERRUPT  H  OOFTOAREERRO  TO  l# 

♦■•ifi  n»«o 

7»*ui  ui  ]Tit«  in4i 
•••SWITCH  ftC16  TAPES 
f««mR 

IOmHAOO  ITftftE  OUTPVT  t/R£SUlTS/T0PVCR\ITR\<4> 

11 ••SWITCH  OUTPUT  VEROP\ITR\<S> 

i2«*6iwc  «cftor\im<4>  iam  cm 

11- *P£STRWY  TAPE1  TAPE2  IAPEJ  TAPE24 
14- •SWITCH  TAPE ft  TAPE! 

13**SUITCH  IIVP  IAPE2 

1  •••SWITCH  ULL  TAPES 
17**XPIX  /  15 

10»«HASS  ST  ME  OUTPVT  t/ftESULTf /T0PI1X\ITIN<4> 
1f**SUlTCH  OUTPUT  01IOP\|Tft\<4> 

2S*»SIV£  01X0P\ITt\<4>  1A04  Ell 

2f*«HASS  STOKE  TAP£20s/RCIULT|/T0AHR\nR\<1X2> 

22*«C0PY  TAPE 20  TP2t\ITR\<4> 

23»*GIVE  ?P20\I1R\<4>  1044  CM» 

24-OESTftlY  TAPE4 
25**SW1TCW  TAPE2  OUT 
24a*$VlTCH  UU17  IAPC4 
27-«XC70  /  2 

28»*MASS  STORE  OUTPWT|/RESW.Tl/TOPC7\fTR\<4> 
2f«*HASS  STORE  TAPE3Sl/R£6WLTSnOftft\ITR\<1X2><3> 
SO»*SUITCH  OUTPUT  C?0P\ITI\<4> 

11»*6IW£  C70PMTRN<4>  1AM  COO 

12- 4SUITCM  TAPES  REIS 
M-*SUITCH  TAPE3  LOLL 

S4>«SUITCM  TAPE4  0U37 
3 •••DESTROY  TAPE!  TAPE2  TAPE10 
••••SWITCH  FI1HT  TAPE1 
1?«»SW1TCH  TAPE 30  TAPE2 
38** SUIT CM  TCRHT  TAPES 
IfMlCII  /  2 

40MRASS  STORE  OUTPUT i/tCtULTS/TOPCOMTR\<4> 

4 1  ••HASS  STORE  TAPEI0t/RCSULTS/TM20MT0\<l><2Xl> 
42**SWITCH  OUTPUT  C00PUTR\<4> 

43**S1VE  COOPMTR\<4>  1AM  COO 
44-»0EST00Y  TAPC1  TAPC2  TAPC1  TAPES 
4S*«SUITCH  TAPES  TCRHT 
4 •■•SWITCH  TCR31  TAPE! 

47MM1TCH  SP<4>  TAPES 


Figure  2.  COSMOS  procedure  file  for  the  experiment 
with  three  iterations 


4I-*SU11CM  TAPEM  Hftl 
49>«tunCH  lit*  TAPES 
5(»XCU  /  4 

S1MNISS  tlMC  OUTPUT  t/RESULTS/T0PCf\ITI\<4> 
S2*»SMITCM  OUTPUT  C90P\ITR\<4> 

53>*I1VE  C90PMTR\<4>  1414  C«» 

94>*SM1TCN  IAPE4  0P<4> 

55- •SWITCH  TAPE!  TCR3T 

94«M*SS  STORE  TAPC2l/*ES«.TS/T0UTC9\ITR\<IX2XI> 

S?-*l£STI0T  1**15  TAPCU  TAPEM  TAPE  SO  TAPEM 

MmUIHN  MV*  TAPES 

S«»*SU!TCH  WK  TAPEM 

4»«*SU1TCH  TAPE2  well 

4M*SUITCH  E1F  TAPE 10 

42**XSMN  /  4 

43»*H*SS  STORE  TAPCMi/RCS0lTS/TITIU\ITt\<1X2Xl> 
44MNASS  ITWE  TAPC4t/RCSUlT|/T9UT0N\IIR\<IX2XI> 
4S**SUITCM  TAPE*  T*\ITA\<4> 

4i«»H«  T4\ITR\<4>  14i4  I.  EM 
4?**SU1TCM  T4\ITR\<4>  TAPE* 

40*»SU!TCN  TAPES  MEN 
4f"*SUITCN  TAPCSO  Elf 
?A**tUlTCH  1API14  Ml 

71MRA8S  STMC  TAPCMs/RCS0LTS/T9PS\ITR\<1X2XI> 
72MI4ITCN  VAPC19  MIT 

73-*KSTRIV  TAPC1  TAPC2  TAPC1I  OUTPUT  TAPC2A  TAPE] 
74*«IF  ITI  .ft.  1  TMti  M  II  U 
?i»* SWITCH  Mil  TAPE* 

74MIC11P  /  2 

77-4MASS  ITMC  OUTPUT  1/OESOLT!/TC1TOP\1TR\<4> 
7S**SU1TCM  OUTPUT  CM9P\!TRV<4> 

7f«*AASS  STMC  TAP£2l/RE8ULTI/TllTF9\lTR\<4> 

M**S«1TCN  IAPC1  TTM 

IlMllE  CT10P\1TR\<4)  1494  CM 

02*'BESTR0T  IAPC1I  MFC 20 

•SmMITCN  TAPC2  FSTAT 

I4>*SUITCM  IISES  TAPE10 

S3-* SUITE*  MAT  IAPE20 

I4«»XC12  /  4 

S7**MASS  SURE  OUTPUT /RE8UITB/TC120P\ITA\<4> 
II>*SU1TCN  OUTPUT  E12QP\1TR\<4> 

IP*»61VE  Ct20P\11R\<4>  MS*  CRB 

90MSS1TCN  TAPE 20  CAM 
fl  ••WITCH  TAPE1S  OHIO 
02*»KtTR01  TAPE  I  TAPES 
93MLCT  ITRP*  ETC  ♦  1 

94MAISS  STORE  TAPC?i/RE6«LTS/IE02CRI\ITRP\<4> 
9SMSUITCN  TAPE?  TAPE  I 


Figure  2.  COSMOS  procedure  file  for  the  experiment 
with  three  iterations  (continued) 


?4«*S0  TO  11 

♦?»*L3*  S8ITCM  HR  TAPE2 
NmKSTIOI  V4KI 

t9«*s«iiCM  inn  urn 

1HMIES1IOT  1APE7  IAPE1 
101**XFCT  /  • 

102**B£STR0T  TAPE8 
103**SUITCH  TAPE2  8862 

104«*HASS  STORE  TAPE7»/R£6ULTS/TFCST<1XIX3> 

105**S8ITCH  TAPE7  FCST<1X2X3> 

104**61V£  FC6T<1X2X3>  1144  R.  CUB 
107**BESTRQY  TAPE! 

108  *<*60  TO  LAS 
1I?»*14*  SO  TO 

110«*LASs  WITCH  FCST<1X2X3>  TAPE1 
111*? 

1I2**HASS 

H3«Gn  Ff  TER31  tTERRSSTHC  HUAtTAPCS  816 (TAPE  14  OHM  ElFtTAPEIOI 
114*8B62:TAPE2  F66CtF66E26712 
KS«l£FAUlT  8IR*/MALTUM6 

114>Gn  CNATsCAATRIXN  REI6xRE1631X12  U8PlU0P8RFIM32  LALL t LALL8UH328 
1 1 7*8837 t 8VPNMF 37LAT32  TERM! iTEKMT 144  SP1:SP1712Z  TAPE1 x6SH07l2 
11 •■SET  B1R*/EXEC8TE 

11?*6EI  XB1X  XFCT  XC4R  X8EI  XC7B  XC86  XCfS  XC11P  XC12  X6NH  XASP  XPTTsXCZ 

120*ENB 

121MXPTT 

122»* SWITCH  TAPE2  F61NI 

123**S8I1CM  TAPES  TTEH 

124«*AASS  STORE  FilNT  t /RESULT 6/T1RTF81 

12S**»ESTR0V  TAP £2  TAPES 

124**68 ITCH  F6GE  TAPE2 

127-*F1L£  NAH£*]WMn 

128*  IOTA  ]BATE*7,ITINE*12  SERB 

12?**S81TCM  TER31  TAPES 

130*»XASP  /  10 

111«*AASS  STORE  TAPEAAl/RESULTS/TASAPI 
132**SUITCH  TAPE!  6SFCT 
133*«DESTR0Y  TAPE10 
134>*S6ITCH  TAPES  TER31 
13S**S81TCN  TAPE44  TAPE1 
134MW1TCM  CHAT  TAPE  II 
137-«XC4R  /  2 

1 38**HASS  STORE  TAPEISs/RCSULTS/TQBSESI 

13?**IESTR0T  TAPE1  TAPE 20 

140**MITCH  TAPE10  CHAT 

141«*SU1TCH  TAPE IS  OISEO 

142**HASS  STORE  TAPESt/RESULTS/TCOZEROU 

143MS81TCH  TAPES  TAPE1 

144**COSHOS  I-C8SITR  8ITH 

Figure  2.  COSMOS  procedure  file  for  the  experiment  run 
with  three  iterations  (continued) 
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